一种区域搜索获取平面光学元件光圈数的方法

提出一种获取平面光学元件光圈数的方法。通过确定一幅干涉条纹图像中三块测试区域及相应的搜索区域,寻找条纹弯曲量的平均值和相应条纹间距的平均值,根据测试区域对应的权重系数,得到被测平面光学元件光圈数的最佳估计值。通过测量标准光圈并比较目视判读结果,得到最大相对误差为5%。实验结果表明,该方法提高了处理干涉条纹图像的适应性,并可替代人工目视判读干涉条纹图像。

MANDIBULAR MORPHOMETRICS AMONG MACAQUES - THE CASE OF MACACCA-THIBETANA

Recent field studies suggest that Macaca thibetana, a large endemic Chinese macaque, may be quite folivorous, distinguishing it from most other macaque species, which tend to be primarily frugivorous. To understand how this diet affects its masticatory system, we conducted a comparative morphometric study of mandibular dimensions. We took linear measurements from male and female mandibles of this species as well as four other macaques-M. fascicularis, M. nemestrina, M. arctoides, and M. assamensis-and four species of Presbytis-P. obscura, P. rubicunda, P. cristata, and P. phayrei-and subjected to them to a variety of analyses. Based on analyses of variances and discriminant analyses on each sex individually, the mandible of M. thibetana corresponds to expected patterns for folivorous primates with respect to its wide condyles and thick corpora: However, the height of the corpus and symphysis are lower, and the anteroposterior length of the condyle is longer than predicted for a folivore. In addition to interpretations specifically relating to M. thibetana, we also discuss the functional morphology of the other species in light of what is published about their diets.

Full counting statistics of transport through two-channel Coulomb blockade systems

A mesoscopic Coulomb blockade system with two transport channels is studied in terms of full counting statistics. It is found that the shot noise and skewness are crucially affected by the quantum mechanical interference. In particular, the super-Poisson behavior can be induced as a consequence of constructive interference, and can be understood by the formation of effective fast-and-slow transport channels. Dephasing and finite temperature effects are carried out together with physical interpretations.

AlGaInAs narrow stripe selective growth on substrates patterned with different mask designs

We have performed a narrow stripe selective growth of oxide-free A1GaInAs waveguides on InP substrates patterned with pairs of SiO2 mask stripes under optimized growth conditions. The mask stripe width varied from 0 to 40 mu m, while the window region width between a pair of mask stripes was fixed at 1.5, 2.5 or 3.5 mu m. Flat and smooth A1GaInAs waveguides covered by specific InP layers are successfully grown on substrates patterned with different mask designs. The thickness enhancement ratio and the photoluminescence (PL) spectrum of the A1GaInAs narrow stripe waveguides are strongly dependent on the mask stripe width and the window region width. In particular, a large PL wavelength shift of 79 nm and a PL FWHM of less than 64 meV are obtained simultaneously with a small mask stripe width varying from 0 to 40 mu m when the window region width is 1.5 mu m. We present some possible interpretations of the experimental observations in considering both the migration effect from a masked region and the lateral vapour diffusion effect.

Size and shape evolution of self-assembled coherent InAs/GaAs quantum dots influenced by seed layer

The size and shape Evolution of self-assembled InAs quantum dots (QDs) influenced by 2.0-ML InAs seed layer has been systematically investigated for 2.0, 2.5, and 2.9-ML deposition on GaAs(1 0 0) substrate. Based on comparisons with the evolution of InAs islands on single layer samples at late growth stage, the bimodal size distribution of InAs islands at 2.5-ML InAs coverage and the formation of larger InAs quantum dots at 2.9-ML deposition have been observed on the second InAs layer. The further cross-sectional transmission electron microscopy measurement indicates the larger InAs QDs: at 2.9-ML deposition on the second layer are free of dislocation. In addition, the interpretations for the size and shape evolution of InAs/GaAs QDs on the second layer will be presented. (C) 2001 Elsevier Science B.V. All lights reserved.

Investigations on oxide-free InGaAlAs waveguides grown by narrow stripe selective MOVPE

Oxide-free InGaAlAs waveguides have been grown on the InP substrates patterned with pairs of SiO2 mask stripes using narrow stripe selective MOVPE. The mask stripe width is varied from 0 to 40 pm, while the window region width between a pair of mask stripes is fixed at 1.5, 2.5 and 3.5 mu m, respectively. Smooth surface s and flat interfaces are obtained in the selectively grown InQaAlAs waveguides. There exhibit strong dependences of the thickness enhancement ratio and the photoluminescence (PL) spectrum on the mask stripe width and the window region width for the InGaAlAs wavegwdes. A large PL peak wavelength shift of 79 nm and a PL full width of at half maximum (FWHM) of less than 64 meV are obtained simultaneously. Some possible interpretations for our investigations are presented by considering both the migration effect from a masked region (MMR) and the lateral vapor diffusion effect (LVD).

Size and shape evolution of self-assembled coherent InAs/GaAs quantum dots influenced by seed layer

The size and shape Evolution of self-assembled InAs quantum dots (QDs) influenced by 2.0-ML InAs seed layer has been systematically investigated for 2.0, 2.5, and 2.9-ML deposition on GaAs(1 0 0) substrate. Based on comparisons with the evolution of InAs islands on single layer samples at late growth stage, the bimodal size distribution of InAs islands at 2.5-ML InAs coverage and the formation of larger InAs quantum dots at 2.9-ML deposition have been observed on the second InAs layer. The further cross-sectional transmission electron microscopy measurement indicates the larger InAs QDs: at 2.9-ML deposition on the second layer are free of dislocation. In addition, the interpretations for the size and shape evolution of InAs/GaAs QDs on the second layer will be presented. (C) 2001 Elsevier Science B.V. All lights reserved.

Sputtering induced by highly charged Pb-208(q+) bombardment on Al surface

Highly charged ions (HCls) carrying high Coulomb potential energy (E-p) could cause great changes in the physical and chemical properties of material surface when they bombard on the solid surface. In our work, the secondary ion yield dependence on highly charged Pbq+ (q = 4-36) bombardment on Al surface has been investigated. Aluminum films (99.99%) covered with a natural oxide film was chosen as our target and the kinetic energy (E-k) was varied between 80 keV and 400 keV. The yield with different incident angles could be described well by the equation developed by us. The equation consists of two parts due to the kinetic sputtering and potential sputtering. The physical interpretations of the coefficients in the said equation are discussed. Also the results on the kinetic sputtering produced by the nuclear energy loss on target Surface are presented.

Seismic profiles across the middle Tan-Lu fault zone in Laizhou Bay, Bohai Sea, eastern China

The stratigraphic architecture, structure and Cenozoic tectonic evolution of the Tan-Lu fault zone in Laizhou Bay, eastern China, are analyzed based on interpretations of 31 new 2D seismic lines across Laizhou Bay. Cenozoic strata in the study area are divided into two layers separated by a prominent and widespread unconformity. The upper sedimentary layer is made up of Neogene and Quaternary fluvial and marine sediments, while the lower layer consists of Paleogene lacustrine and fluvial facies. In terms of tectonics, the sediments beneath the unconformity can be divided into four main structural units: the west depression, central uplift, east depression and Ludong uplift. The two branches of the middle Tan-Lu fault zone differ in their geometry and offset: the east branch fault is a steeply dipping S-shaped strike-slip fault that cuts acoustic basement at depths greater than 8 km, whereas the west branch fault is a relatively shallow normal fault. The Tan-Lu fault zone is the key fault in the study area, having controlled its Cenozoic evolution. Based on balanced cross-sections constructed along transverse seismic line 99.8 and longitudinal seismic line 699.0, the Cenozoic evolution of the middle Tan-Lu fault zone is divided into three stages: Paleocene-Eocene transtension, Oligocene-Early Miocene transpression and Middle Miocene to present-day stable subsidence. The reasons for the contrasting tectonic features of the two branch faults and the timing of the change from transtension to transpression are discussed. Crown Copyright (C) 2008 Published by Elsevier Ltd. All rights reserved.

Distribution characteristics of matrix-bound phosphine along the coast of China and possible environmental controls

Matrix-bound phosphine (MBP) concentrations in surface sediments collected from 37 stations along the coast of China in 2006 are reported. MBP was found in all samples and the average concentration was 6.30 ng kg(-1) dry weight (dw). The distribution of MBP showed certain spatial variation characteristics with high MBP concentrations at stations near to the coast. The average concentrations of MBP in the northern Yellow Sea (NYS), the southern Yellow Sea (SYS), the northern area of East China Sea (NECS), the southern area of East China Sea (SECS), and South China Sea (SCS) were 5.57 +/- 3.78, 3.78 +/- 2.81, 5.27 +/- 3.07, 5.48 +/- 4.05 and 13.52 +/- 7.86 ng kg(-1) dw. respectively. The correlations between MBP and influencing factors, such as the sedimentary environmental characteristics (sediment type, the grain size, contents of phosphorous, organic matters and redox potential) and the aquatic environmental characteristics (temperature, salinity, depth and hydrodynamics) were studied. The results indicated that MBP was strongly influenced by various factors, such as total phosphorus (TP), organic phosphorus (OP), organic carbon (OC), the grain size and hydrodynamics, all of which not only offered reasonable interpretations for the distribution characteristics of MBP but also provided evidence to support the viewpoint that phosphine originated from OP decomposition. This work is the first comprehensive study of the distribution of MBP along the coast of China and its relationships with environmental factors which will lead to a better understanding of the phosphorus (P) biogeochemical cycle in the sea. (C) 2008 Elsevier Ltd. All rights reserved.

地震波传播模拟的边界元―体积元积分方程法及其近似求解

To improve the efficiency of boundary-volume integral equation technique, this paper is involved in the approximate solutions of boundary-volume integral equation technique. Firstly, based on different interpretations of the self-interaction and extrapolation operators of the resulting boundary integral equation matrix, two different hybrid BEM+Born series modeling schemes are formulated and validated through comparisons with the full-waveform BE numerical solutions for wave propagation simulation in a semicircular alluvial valley and a complex fault model respectively. Numerical experiments indicate that both the BEM+Born series modeling schemes are suitable for complex geological structures and significantly improve computational efficiency especially for the cases of high frequencies and multisource seismic survey. Then boundary-volume integral equation technique is illuminated in detail and verified by modeling wave propagation in complex media. Furthermore, the first-order and second-order Born approximate solutions for the volume-scattering waves are studied and quantified by numerical simulation in different random medium models. Finally, preconditioning generalized minimal residual method is applied to solve boundary-volume integral equation and compared with Gaussian elimination method. Numerical experiments indicate this method makes the calculations more efficient.

东海陆架盆地西南部烟囱构造与油气运聚关系研究

In recent years, chimney structure has been proved one of important indicators and a useful guide to major petroleum fields exploration through their exploration history both at home and abroad. Chimney structure, which has been called "gas chimney" or "seismic chimney", is the special fluid-filled fracture swarm, which results from the boiling of active thermal fluid caused by abruptly decreasing of high pressure and high temperature in sedimentary layers of upper lithosphere. Chimney structure is well developed in continental shelf basin of East China Sea, which indicates the great perspectives of petroleum resources there. However, the chimney structure also complicated the petroleum accumulation. So the study of chimney structure on its formation, its effect on occurrence and distribution of petroleum fields is very important not only on theoretical, but also on its applied research. It is for the first time to make a clear definition of chimney structure in this paper, and the existence and practical meaning of chimney structure are illustrated. Firstly, on the viewpoint of exploration, this will amplify exploration area or field, not only in marine, but also on continent. Secondly, this is very important to step-by-step exploration and development of petroleum fields with overpressure. Thirdly, this will provide reference for the study on complex petroleum system with multi-sources, commingled sources and accumulation, multi-stage accumulations, and multi-suits petroleum system in the overlay basin. Fourthly, when the thermal fluid enters the oceanic shallow layer, it can help form gas hydrate under favorable low-temperature and high-pressure conditions. Meanwhile, the thermal fluid with its particular component and thermal content will affect the physical, chemical and ecological environments, which will help solving the problem of global resources and environment. Beginning from the regional tectonic evolution characteristics, this paper discussed the tectonic evolution history of the Taibei depression, then made an dynamical analysis of the tectonic-sedimentary evolution during the Mesozoic and Cenozoic for the East China Sea basin. A numerical model of the tectonic-thermal evolution of the basin via the Basin-Mod technique was carried out and the subsidence-buried history and thermal history of the Taibei depression were inverse calculated: it had undergone a early rapid rift and sag, then three times of uplift and erosion, and finally depressed and been buried. The Taibei depression contains a huge thick clastic sedimentary rock of marine facies, transitional facies and continental facies on the complex basement of ante-Jurassic. It is a part of the back-arc rifting basins occurred during the Mesozoic and Cenozoic. The author analyzed the diagenesis and thermal fluid evolution of this area via the observation of cathodoluminescence, scanning electron microscope and thin section, taking advantage of the evidences of magma activities, paleo-geothermics and structural movement, the author concluded that there were at least three tectonic-thermal events and three epochs of thermal-fluid activities; and the three epochs of thermal-fluid activities were directly relative to the first two tectonic-thermal events and were controlled by the generation and expulsion of hydrocarbon in the source rock simultaneously. Based on these, this paper established the corresponding model between the tectonic-thermal events and the thermal-fluid evolution of the Taibei Depression, which becomes the base for the study on the chimney structures. According to the analyses of the gas-isotope, LAM spectrum component of fluid inclusion, geneses of CO_2 components and geneses of hydrocarbon gases, the author preliminarily verified four sources of the thermal fluid in the Taibei Depression: ① dehydration of mud shale compaction, ② expulsion of hydrocarbon in the source rock; ③ CO_2 gas hydro-thermal decomposition of carbonatite; ④magma-derived thermal fluid including the mantle magma water and volatile components (such as H_2O, CO_2, H_2S, SO_2, N_2 and He etc.). On the basis of the vitrinite reflectance (Ro), homogenization temperature of fluid inclusion, interval transit time of major well-logging, mud density of the wells, measured pressure data and the results of previous studies, this paper analyzed the characteristics of the geothermal fields and geo-pressure fields for the various parts in this area, and discussed the transversal distribution of fluid pressure. The Taibei depression on the whole underwent a temperature-loss process from hot basin to cold basin; and locally high thermal anomalies occurred on the regional background of moderate thermal structure. The seal was primarily formed during the middle and late Paleocene. The overpressured system was formed during the middle and late Eocene. The formation of overpressured system in Lishui Sag underwent such an evolutionary process as "form-weaken-strengthen-weaken". Namely, it was formed during the middle and late Eocene, then was weakened in the Oligocene, even partly broken, then strengthened after the Miocene, and finally weakened. The existence of the thermal fluid rich in volatile gas is a physical foundation for the boiling of the fluid, and sharply pressure depletion was the major cause for the boiling of the fluid, which suggests that there exists the condition for thermal fluid to boil. According to the results of the photoelastic simulation and similarity physical experiments, the geological condition and the formation mechanism of chimnestructures are summarized: well compartment is the prerequisite for chimney formation; the boiling of active thermal fluid is the original physical condition for chimney formation; The local place with low stress by tension fault is easy for chimney formation; The way that thermal fluid migrates is one of the important factors which control the types of chimney structures. Based on where the thermal fluid come from and geometrical characteristics of the chimney structures, this paper classified the genetic types of chimney structures, and concluded that there existed three types and six subtypes chimney structures: organic chimney structures generated by the hydrocarbon-bearing thermal fluid in middle-shallow layers, inorganic and commingling-genetic chimney structures generated by thermal fluid in middle-deep layers. According to the seismic profiles interpretations, well logging response analysis and mineralogical and petrological characteristics in the study area, the author summarized the comprehensive identification marks for chimney structures. Especially the horizon velocity analysis method that is established in this paper and takes advantage of interval velocity anomaly is a semi-quantitative and reliable method of chimney structure s identification. It was pointed out in this paper that the occurrence of the chimney structures in the Taibei depression made the mechanism of accumulation complicated. The author provided proof of episodic accumulation of hydrocarbon in this area: The organic component in the boiling inclusion is the trail of petroleum migration, showing the causality between the boiling of thermal fluid and the chimney structures, meanwhile showing the paroxysmal accumulation is an important petroleum accumulation model. Based on the evolutionary characteristics of various types of chimney structures, this paper discussed their relationships with the migration-accumulation of petroleum respectively. At the same time, the author summarized the accumulating-dynamical models associated with chimney structures. The author analyzed such accumulation mechanisms as the facies state, direction, power of petroleum migration, the conditions of trap, the accumulation, leakage and reservation of petroleum, and the distribution rule of petroleum. The author also provides explanation for such practical problems the existence of a lot of mantle-derived CO_2, and its heterogeneous distribution on plane. By study on and recognition for chimney structure, the existence and distribution of much mantle-derived CO_2 found in this area are explained. Caused by tectonic thermal activities, the deep magma with much CO_2-bearing thermal fluid migrate upward along deep fault and chimney structures, which makes two wells within relatively short distance different gas composition, such as in well LF-1 and well LS36-1-1. Meanwhile, the author predicted the distribution of petroleum accumulation belt in middle-shallow layer for this area, pointed out the three favorable exploration areas in future, and provided the scientific and deciding references for future study on the commingling-genetic accumulation of petroleum in middle-deep layer and the new energy-gas hydrate.

各向异性介质中的地震面波和S波分裂

Surface wave propagation in the anisotropic media and S-wave splitting in China mainland are focused in this M.S. dissertation. We firstly introduced Anderson parameters in the research of surface wave propagation in the anisotropic media were deduced, respectively. By applying the given initial model to the forward calculation of Love wave, we compared dispersion curves of Love wave in the anisotropic media with the one in the isotropic media. the results show that, although the two kind of results are similar with each other, the effect of anisotropy can not be neglected. Furthermore, the variation of anisotropy factors will result in the variation of dispersion curves, especially for high-mode one. The method of grid dispersion inversion was then described for further tectonic inversion. We also deduced inversion equation on the condition that the layered media is anisotropic, and calculated the phase-velocity partial derivatives with respect to the model parameters, P- and S-wave velocities, density, anisotropic parameters for Rayleigh wave and Love wave. Having analyzed the results of phase-velocity partial derivatives, we concluded that the derivatives within each period decreased with the depth increasing, the phase-velocity of surface wave is sensitive to the S-wave velocities and anisotropic factors and is not sensitive to the densities of layers. Dispersion data of Love wave from the events occurred during the period from 1991 to 1998 around the Qinghai and Tibet Plateau, which magnitudes are more than 5.5, have been used in the grid dispersion inversion. Those data have been preprocessed and analyzed in the F-T domain. Then the results of 1°*1° grid dispersion inversion, the pure path dispersion data, in the area of Qianghai and Tibet Plateau were obtained. As an example, dispersion data have been input for the tectonic inversion in the anisotropic media, and the results of anisotropic factors under the region of Qianghai and Tibet Plateau were initially discussed. As for the other part of this dissertation. We first introduced the phenomena of S-wave splitting and the methods for calculation the splitting parameters. Then, We applied Butterworth band-pass filter to S-wave data recorded at 8 stations in China mainland, and analyzed S-wave splitting at different frequency bands. The results show the delay time and the fast polarization directions of S-wave splitting depend upon the frequency bands. There is an absence of S-wave splitting at the station of Wulumuqi (WMQ) for the band of 0.1-0.2Hz. With the frequency band broaden, the delay time of S-wave splitting decreases at the stations of Beijing (BJI), Enshi (ENH), Kunming (KMI) and Mudanjiang (MDJ); the fast polarization direction at Enshi (ENH) changes from westward to eastward, and eastward to westward at Hailaer (HIA). The variations of delay time with bands at Lanzhou (LZH) and qiongzhong (QIZ) are similar, and there is a coherent trend of fast polarization directions at BJI, KMI and MDJ respectively. Initial interpretations to the results of frequency band-dependence of S-wave splitting were also presented.

港西地区河流相砂岩油藏剩余油分布研究

The Gangxi oil field has reached a stage of high water production. The reservoir parameters, such as reservoir physical characteristics, pore structure, fluid, have obviously changed. This thesis therefore carries out a study of these parameters that control reservoir characteristics, physical and chemical actions that have taken place within the reservoirs due to fluid injection, subsequent variations of reservoir macroscopic physical features, microscopic pore structures, seepages, and formation fluid properties. This study rebuilds a geologic model for this oil field, establishes a log-interpreting model, proposes a methodology for dealing with large pore channels and remnant oil distribution, and offers a basis for effective excavation of potential oil, recovery planning, and improvement of water-injection techniques. To resolve some concurrent key problems in the process of exploration of the Gangxi area, this thesis carries out a multidisciplinary research into reservoir geology, physical geography, reservoir engineering, and oil-water well testing. Taking sandstone and flow unit as objects, this study establishes a fine geologic model by a quantificational or semi-quantificational approach in order to understand the remnant oil distribution and the reservoir potential, and accordingly proposes a plan for further exploration. By rebuilding a geological model and applying reservoir-engineering methods, such as numerical simulation, this thesis studies the oil-water movement patterns and remnant-oil distribution, and further advances a deployment plan for the necessary adjustments and increase of recoverable reserves. Main achievements of this study are as follows: 1. The Minghazhen Formation in the Gangxi area is featured by medium-sinuosity river deposits, manifesting themselves as a transitional type between typical meandering and braided rivers. The main microfacies are products of main and branch channels, levee, inter-channel overflows and crevasse-splay floodplains. The Guantao Group is dominantly braided river deposit, and microfacies are mainly formed in channel bar, braided channel and overbank. Main lithofacies include conglomerate, sandstone, siltstone and shale, with sandstone facies being the principal type of the reservoir. 2. The reservoir flow unit of the Gangxi area can be divided into three types: Type I is a high-quality heterogeneous seepage unit, mainly distributed in main channel; Type II is a moderate-quality semi-heterogeneous seepage unit, mainly distributed in both main and branch channels, and partly seen within inter-channel overflow microfacies; Type III is a low-quality, relatively strong heterogeneous seepage unit, mainly distributed in inter-channel overflow microfacies and channel flanks. 3. Flow units and sedimentary microfacies have exerted relatively strong controls on the flowing of underground oil-water: (1) injection-production is often effective in the float units of Type I and II, whilst in the same group of injection-production wells, impellent velocity depends on flow unit types and injection-production spacing; (2) The injection-production of Type III flow unit between the injection-production wells of Type I and II flow units, however, are little effective; (3) there can form a seepage shield in composite channels between channels, leading to inefficient injection and production. 4. Mainly types of large-scale remnant-oil distribution are as follows: (1) remnant oil reservoir of Type III flow unit; (2) injection-production well group of remnant oil area of Type III flow unit; (3) remnant oil reservoirs that cannot be controlled by well network, including reservoir featured by injection without production, reservoir characterized by production without injection, and oil reservoir at which no well can arrive; (4) remnant oil area where injection-production system is not complete. 5. Utilizing different methods to deal with different sedimentary types, sub-dividing the columns of up to 900 wells into 76 chronostratigraphic units. Four transitional sandstone types are recognized, and contrast modes of different sandstone facies are summarized Analyzing in details the reservoirs of different quality by deciphering densely spaced well patterns, dividing microscopic facies and flow units, analyzing remnant oil distribution and its effect on injection-production pattern, and the heterogeneity. Theory foundation is therefore provided for further excavation of remnant oil. Re-evaluating well-log data. The understanding of water-flood layers and conductive formations in the Gangxi area have been considerably improved, and the original interpretations of 233 wells have changed by means of double checking. Variations of the reservoirs and the fluid and formation pressures after water injection are analyzed and summarized Studies are carried out of close elements of the reservoirs, fine reservoir types, oil-water distribution patterns, as well as factors controlling oil-gas enrichment. A static geological model and a prediction model of important tracts are established. Remaining recoverable reserves are calculated of all the oil wells and oil-sandstones. It is proposed that injection-production patterns of 348 oil-sandstones should be adjusted according to the analysis of adaptability of all kinds of sandstones in the injection-production wells.

基于MRF模型的油气储层属性随机建模方法研究

Stochastic reservoir modeling is a technique used in reservoir describing. Through this technique, multiple data sources with different scales can be integrated into the reservoir model and its uncertainty can be conveyed to researchers and supervisors. Stochastic reservoir modeling, for its digital models, its changeable scales, its honoring known information and data and its conveying uncertainty in models, provides a mathematical framework or platform for researchers to integrate multiple data sources and information with different scales into their prediction models. As a fresher method, stochastic reservoir modeling is on the upswing. Based on related works, this paper, starting with Markov property in reservoir, illustrates how to constitute spatial models for catalogued variables and continuum variables by use of Markov random fields. In order to explore reservoir properties, researchers should study the properties of rocks embedded in reservoirs. Apart from methods used in laboratories, geophysical means and subsequent interpretations may be the main sources for information and data used in petroleum exploration and exploitation. How to build a model for flow simulations based on incomplete information is to predict the spatial distributions of different reservoir variables. Considering data source, digital extent and methods, reservoir modeling can be catalogued into four sorts: reservoir sedimentology based method, reservoir seismic prediction, kriging and stochastic reservoir modeling. The application of Markov chain models in the analogue of sedimentary strata is introduced in the third of the paper. The concept of Markov chain model, N-step transition probability matrix, stationary distribution, the estimation of transition probability matrix, the testing of Markov property, 2 means for organizing sections-method based on equal intervals and based on rock facies, embedded Markov matrix, semi-Markov chain model, hidden Markov chain model, etc, are presented in this part. Based on 1-D Markov chain model, conditional 1-D Markov chain model is discussed in the fourth part. By extending 1-D Markov chain model to 2-D, 3-D situations, conditional 2-D, 3-D Markov chain models are presented. This part also discusses the estimation of vertical transition probability, lateral transition probability and the initialization of the top boundary. Corresponding digital models are used to specify, or testify related discussions. The fifth part, based on the fourth part and the application of MRF in image analysis, discusses MRF based method to simulate the spatial distribution of catalogued reservoir variables. In the part, the probability of a special catalogued variable mass, the definition of energy function for catalogued variable mass as a Markov random field, Strauss model, estimation of components in energy function are presented. Corresponding digital models are used to specify, or testify, related discussions. As for the simulation of the spatial distribution of continuum reservoir variables, the sixth part mainly explores 2 methods. The first is pure GMRF based method. Related contents include GMRF model and its neighborhood, parameters estimation, and MCMC iteration method. A digital example illustrates the corresponding method. The second is two-stage models method. Based on the results of catalogued variables distribution simulation, this method, taking GMRF as the prior distribution for continuum variables, taking the relationship between catalogued variables such as rock facies, continuum variables such as porosity, permeability, fluid saturation, can bring a series of stochastic images for the spatial distribution of continuum variables. Integrating multiple data sources into the reservoir model is one of the merits of stochastic reservoir modeling. After discussing how to model spatial distributions of catalogued reservoir variables, continuum reservoir variables, the paper explores how to combine conceptual depositional models, well logs, cores, seismic attributes production history.